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Abstract 

Non-axisymmetric accretion discs can lead to flux variability. Here we provide an 
analytic framework for modeling non-axisymmetric structures like hotspots and spi- 
ral waves and their influence on observed timing measurements. The presently un- 
explained Low- frequency Quasi-Periodic Oscillations (LFQPO), observed in X-ray 
binaries and cataclysmic variables, could be the result of such discs. Our framework 
serves as a guide to quantify the properties that non-axisymmetric structures pro- 
duced by nonlinear accretion disc models must have in order to explain observed 
features such as LFQPOs. The results from our microquasar appliations also provide 
analogous predictions for X-ray modulation in active galactic nuclei. The formalism 
and physical interpration is of practical use for generic non-axisymmetric accretion 
disc systems. 
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1 Introduction 



X-ray emission from accreting black holes in binary stellar systems varies on 
time scales ranging from milliseconds to years. Variability on time scales longer 
than days appears to be driven by changes in the accretion rate onto the black 
hole, and is often manifested as transient outbursts in which the luminosity of 
a source changes by a million- fold Lasota (2001). At the shortest time scales, 
quasi-periodic oscillations (QPOs) are observed in the X-ray emission. The 
highest frequency QPOs (> 100Hz) are consistent with those expected from 
general relativistic orbits near the innermost stable orbit around the black 
hole, and they are likely caused by inhomogeneities in the inner accretion 
flow (Remillard et al., 2002, Stella & Vietri, 1998). However, the cause of 
low frequency QPO (LFQPO)s 0.1-20 Hz is still a mystery. This motivates 
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the work herein. For neutron star X-ray binaries, low frequency QPOs can 
be sub-divided into more precise categories (e.g. Psaltis et al.,1999). Here 
we focus on a basic paradigm for the modulation and leave explanations of 
possible harmonic relations and phase lag behavior Varniere (2005) for further 
work. Black hole microquasars avoid the role of any solid stellar surface, and 
offer a purer probe of the accretion process without interaction with a stellar 
surface. The basic properties of LFQPOs that any model needs to explain are 
these: 

(1) LFQPOs often appear to be present in widely spaced observations, over 
a period of weeks, with fractionally narrow frequency widths (Au/u ~ 1/30) 
(Morgan et al., 1997, Remillard et al., 1999). Although the observations are 
not continuous, the fact that the LFQPO are observed to be the same (within 
a very small variation) over several observations within the same week/month, 
suggest that the QPO are likely always present during that time, at least when 
averaged over such long time scales within the hard state; there are periods 
where the LFQPO appear and disappear on timescale of a few seconds. If 
LFQPOs result from inhomogeneities orbiting at a Keplerian speed in the 
accretion disc they are constrained to a narrow range of radii. 

(2) LFQPO amplitudes are typically 5-10% RMS (root mean squared), but 
can reach 20% RMS. Given that the X-ray luminosities of black hole binaries 
approach 10 39 erg s" 1 , LFQPOs involve an enormously energetic fraction of 
the accretion flow. 

(3) On the other hand, LFQPOs are transient features, so the spectral proper- 
ties of the X-ray emission when LFQPOs are present and absent can be used 
to constrain the LFPQO origin (Muno et al, 1999, Sobczak et al, 2000). 

(4) The 0.1-20 Hz LFQPOs only appear when the non-thermal component 
of the X-ray spectrum is strong, and their fractional amplitudes increase with 
energy between 2-20 keV. However, a thermal component that contributes 
~ 10% of the X-ray emission at lower energies must also be present, and the 
frequencies of the LFQPO appear to be correlated best with this. The energy 
output of accreting black holes can generally be decomposed into a ~ 1 keV 
thermal component (thought to originate from the optically thick accretion 
disc) and a non-thermal component with a power-law spectrum extending 
beyond 100 keV (thought to result from inverse- Compton scattering of cool 
photons by a corona of hot electrons.) Thus, it appears that the 0.1-20 Hz 
QPOs either originate in the boundary between the disc and the corona, or 
are part of the mechanism which accelerates the hot Comptonizing electrons. 

As mentioned, LFQPOs can be further subdivided (Psaltis et al., 1999), but 
the above represents a basic set of characteristics that a zeroth order LFQPO 
flux modulation model should explain. 
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Two different conceptual paradigms for LFQPOs have been proposed in this 
regard: (1) In the centrifugal pressure supported boundary layer model (CEN- 
BOL, Chakrabarti & Manickam, 2000), a QPO is produced by a shock in the 
accretion flow where it makes a transition from a Keplerian disc to a hot 
Comptonizing region. If the cooling time of the post-shock region is resonant 
with the free-fall time at the shock, the shock can oscillate radially with a fre- 
quency on order a Hz. The QPO is, in principle, produced because the shock 
modulates the flux of seed photons that reach the Comptonizing post-shock 
region. In this model, the LFQPO would represent a global radial oscillation 
that modifies the emitted flux. (2) A second paradigm for LFQPOs is non- 
axisymmetric structures in the disk. Their motion around the black hole would 
lead to flux modualtion. These structures could move at the Keplerian speed, a 
precession speed (such as Lense-Thirring) or a phase velocity associated with a 
structural instability. One example occurs in the accretion-ejection instability 
model (AEI, Tagger & Pellat, 1999, Varniere & Tagger, 2002). Here a spiral 
shock forms in an accretion disc threaded by a vertical magnetic field. The dis- 
persion relation is similar to that of a galactic disc but with the gravitational 
potential replaced by a magnetic one. At the co-rotation radius, at which the 
Keplerian velocity matches that of the pattern, a Rossby vortex forms and 
emits vertical Alfven waves into a corona. The LFQPO would result from the 
orbit of this non-axisymmetric structure in the disc. 2D numerical simulations 
of this instability, for example, demonstrated (Varniere et al, 2003) that the 
X-ray flux was modulated. However, in 2D, the constraints on the cause of the 
modulation and how it can be stronger and weaker were limited. 



In this paper, we focus on the paradigm exemplified by the latter type of 
model, namely the production of LFQPOs by the orbit of a non-axisymmetric 
structure. Here we do not focus on the specific origin of the structure but 
rather investigate the consequences of the presence of such a structure. We 
investigate the specific question of whether a stable pattern in an accretion disc 
can reproduce the observed characteristics of LFQPO simply by its motion 
around the central engine. In Sec. 2 we give analytical formulae that can 
be used to model arbitrarily shaped blobs and spirals as non-axisymmetric 
disc features. In Sec. 3 we compute the emission from an accretion disc with 
such non-axisymmetric structures and discuss the influence of the shape of 
the particular non-axisymmetric structures on the RMS LFQPO amplitude. 
We also discuss the parameter choices which provide the most observationally 
consistent non-axisymmetric structures for the microquasar example discussed 
above. We also note that the formalism is generically applicable to other non- 
axisymmetric accretion disc systems. We conclude in Sec. 5. 
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2 Modeling non-axisymmetric discs 

Blobs and spirals are two useful categories of non-axisymmetric structures. We 
use blob to indicate a generic localized feature and spirals to indicate the more 
specifically identifiable global structure produced by gravitational or MHD 
instabilities-the latter via the accretion-ejection instability (AEI) i.e. Tagger 
& Pellat 1999. Here we do not detail the formation of these non-axisymmetric 
structures but focus on providing an analytic framework that characterizes 
their shape for practical use and allows observational implications for flux 
modulation to be quantified. 

Because discs of observed systems such as microquasars and active galactic 
nuclei (AGN) are not resolved, the disc structure cannot be directly imaged. 
But from studying the timing evolution of the flux, especially the presence 
of the LFQPO with its frequency and RMS amplitude, we can infer what 
structures might be present. We concentrate on the RMS amplitude of the 
modulation, assuming that the frequency is already matched by the presumed 
location of the given structure. 

Taking into account the disc thickness is a particularly important aspect of 
our endeavor. For discs viewed at highly inclined angles, the shadowing from 
a local thickening of the disc can be important. The height and temperature 
profiles are coupled in the hydrostatic equilibrium approximation. Although 
we make this approximation here, detailed 3D simulations will ultimately be 
needed to get an exact profile of h and T for a more realistic disc. 

2.1 Analytic expression for the disc thickness 
For the disc thickness, we write 

h(r, 0) = h Q (r) + /*i(r(0), 4>)) = h (r) + s(r - r s )d(r) (1) 

where hi is a perturbation in thickness around the unperturbed thickness h Q , 
r is the radial location, and <fi is the azimuthal angle. We have used h\ = 
s(r — r s )d(r), where s is a "shape" function. The latter is finite only near the 
disc structure causing the non-axisymmetry. The blob or spiral wave feature 
is localized by r s = r c e Q<r ^, where r c is the point where the structure begins 
and a is the opening angle of the structure. The quantity d(r) is a thickness 
function which is defined as the height of the disc at each point. 

We now take s(r — r s ) to be a Gaussian and d to be a power-law in r. This 
provides a simple but useful framework to model non-axisymmetric structures. 
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From (1) we then have 



h(r, 0) = h (r) + 7 



r, 



r 



c 



) 



(2) 



where the constant (3 measures how fast the thickness decreases from the 
maximum, 8 parameterizes the radial extent of the structure, and 7 defines 
the disc thickness at r c , We allow the maximum height 7 to be a function 
of the unperturbed thickness at that point, that is 7 = jh (r c ). This allows 
consideration of cases with similar ho/r but different r c . We also consider the 
number of times the non-axisymmetric structure winds around the disc. 

The influence of the parameters in (2) is illustrated in Fig. 1. The role of a is 
seen in Fig. la which shows a top view of two spirals' "spines" (i.e. the line 
defined by r s (a, 0, r c ) tracing the maximal height above the unperturbed disc 
at each 0.). The larger the a the more open the spiral. The role of parameters 
7, 5, and (3 are seen in Fig. lb which shows cross sections of the disc height 
profile for different projected into the r,z plane for r c = 1.5. Increasing S 
would increase the full width-half-maximum of the peaks, which in practice 
has a larger effect on the rise to the peak than the fall because the Gaussian 
perturbation is superimposed upon a positively sloped disk. Increasing 7 would 
increase the maximum height above the h Q , and the downward slope of the 
line that connects the peaks at different values would be steeper for larger (3. 
Fig. lc shows a 3-D close-up of the inner region of the disc with a one-armed 
spiral characterized by h Q /r = 0.01, r c = 3, a = 0.07, 7 = 0.3, (3 = 2, and 
5 = 0.2 at an viewing angle of 70° from the normal. 



2.2 Emission 



Having obtained an expression for the disc height, we compute the temperature 
using the approximation that c s = hVL, which gives T = /// R gas c 2 . This means, 
for example, that a change in temperature by 20% is related to a change in 
thickness of about 45%. 

In order to determine the observational effect of a non-axisymmetric disc thick- 
ness, we compute the flux as a function of azimuth that an observer would 
received as time-dependent flux modulation when the disc rotates. Assum- 
ing that the spectrum from each point is a blackbody, we use the height- 
temperature relation above to compute the photon flux 



f{E) 



2 E 2 



(3) 



c 2 h 3 exp[E/kT] - 1' 
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We then sum the flux from each cell, multiply by cos([tt/2 — 9] — (), where 9 
is the disc inclination angle from the normal and 



We trace rays from the individual cells to the observer in order to determine 
whether the flux is intercepted by the outer portions of the disc. We compute 
the observed spectrum as the structure rotates by viewing a single snapshot 
at all angles 0. The spectrum is taken as multi-temperature blackbody, and 
we sum the number of photons received. We compute the Fourier transform of 
the profile in order to estimate the amplitudes of any modulation at multiples 
of the pattern frequency. 

Our simple spectral modeling does not include scattering in the disc atmo- 
sphere, the presence of an electron corona, or special and general relativistic 
(GR) effects in part because the LFQPOs observed in microquasars could be 
coming from shadowing by structures at > 100 gravitational radii. Also, the 
non-relativistic formalism applies at all radii for generic disks systems around 
stars. We note however, that non-Keplerian origins of the non-axisymmetries 
leading to LFQPOs are also possible, such as binary induced precession, spiral 
wave phase velocities that differ from Kelperian, or Lense-Thirring precession, 
The latter is an intrinsically general relativistic phenomenon. While relativistic 
effects should be considered in future more detailed applications to compact 
systems, here we focus on qualifying the effect of the key parameters in (2) on 
the modulation, not on the particular source of the the non-axisymmetry. 



3 Simulated flux modulation 



Table 1 shows a subset of our numerically solved cases for the blob and spiral. 
Cases 1-15 can be considered spirals, whilst cases 16-19 can be considered 
blobs because their structures extend less than 2n radians. For all cases we 
have taken an inclination angle of 70°, motivated by the inferred observation 
angle of GRS 1915+105 (Mirabel & Rodrguez, 1994). 

There is degeneracy in the relative influences of the various parameters of (2) 
on the RMS amplitude of modulation of Table 1, but it is instructive to discuss 
the influence of each of the parameters separately. 




(4) 
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3. 1 Spirals and physical interpretation of the geometric parameters 

First notice the lack of influence of the range of > 2tt in Table 1: For each 0, 
most of the modulation comes from the most inward "bump" in h(r). Allowing 
for more than 1 bump does not modify the RMS modulation significantly, as 
seen in cases #1, #2, #3. For the majority of the other cases we will therefore 
consider only the range of <f) < 2n. 

The parameter 7, measures the maximum thickness of the perturbation. Its 
effect is evident in comparing cases #3, #6, #7. The greater 7, all else being 
equal, the stronger the modulation. This is because a thicker structure more 
strongly shadows the inner disc. A change in thickness of 70% means a temper- 
ature change of ~ 50%. This somewhat extreme case is taken to illustrate the 
role of 7 on the RMS amplitude. Constraining the amplitude modulation from 
first principles requires a 3D MHD disc simulation that exhibits a spiral wave 
instability. Here the non-axisymmetry would emerge rather than be imposed. 
Previous 2.5-D (i.e. no dynamical vertical structure) disc simulations unstable 
to the AEI exhibit spiral waves, and lead to a temperature variation typically 
of order of 20 — 30% in the hydrostatic equilibrium approximation (Varniere 
et al, 2003). But in 3-D, hydrostatic equilibrium would underestimate the per- 
turbation thickness of an AEI generated spiral because of additional heating 
from spiral shock dissipation. Thus a higher 7 than that obtained from the 
2.5D simulation can be expected, provided the AEI survives in 3-D. 

Because 7 = 7/i (r c ) is the maximum height of the spiral arm, a change in h 
(the zeroth order disc thickness) also increases the modulation. By comparing 
cases #14 and #15 we see that changing h Q by a factor of two is not exactly 
the same as changing 7 by a factor of two because increasing 7 also increases 
the difference between the maximum height of the perturbation and h Q . 

As discussed at the end of Sec. 2, the parameter f3 measures how fast the 
maximum height along the spiral decreases with radius. Its effect is revealed 
in cases #3, #4, #5. The RMS amplitude of the modulation increases with 
(3 for an a > spiral. A rapidly decreasing height perturbation with radius 
means that the height also rapidly decreases along the spiral. This causes 
a stronger modulation by producing a stronger azimuthal variation. This is 
particularly important for tight spirals (small a ~ 0.05): were it not for a large 
f3, little non-axisymmetry would otherwise arise. 

The simulations relevant for studying the influence of r c , the initial radius 
where the perturbation begins (r c ), are #6, #9, #10 and #11. We see a 
maximum in the RMS amplitude for r c around 2 rj n , where r« n is the inner 
disc radius: When the spiral is too near r in , it obscures less of the inner disc 
which is the most luminous part. On the other hand, if the spiral is too far out, 
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it can only obscure a less luminous outer region giving a smaller modulation. 
The combination of these effects is an intermediate r c for maximal modulation. 
Since the position of the spiral is chosen based on the LFQPO frequency, a 
correlation between the RMS amplitude and the frequency of the QPO is 
expected. In practice this may be hard to disentangle from variations in the 
other parameters. 



The influence of the parameter a, which determines the opening of the spiral 
wave can be seen from simulations #6, #12 and #13. There we Ssee that 
the more open the spiral, the higher the RMS amplitude. This is because a 
more open spiral means a more non-axisymmetric thickness profile, leading to 
a higher RMS amplitude of modulation. 



Finally, consider the parameter 8, which measures the width of the spiral or 
blob at its base. From cases #6 and #8 we see that a larger 5, implies a 
larger RMS amplitude. For our choice of a highly inclined system, this trend 
results because a larger 5 makes more of the region inner to the peak of the 
perturbation more perpendicular to the line of sight. This produces a larger 
observed flux when the observer is looking at the disk from an azimuth for 
which the line of sight intersects the inner part of the perturbation. But as 
the disk rotates, the outer edge of the perturbation comes into view, and 
the shadowing of the inner region occurs similarly for large or small 8. The 
contrast in flux (and thus the RMS amplitude of the modulation) is thus larger 
for larger 8, explaining the trend. 



3.2 Blobs or hot spots 



Because the difference between spirals and blobs in our model is just whether 
the structure extends for a range of > 27r (spiral) or < 2n (blob), the 
parameter influences are similar for the two cases. The relevant blob simula- 
tions are #16, #17, #18, and #19. There we see that a blob can create a 
non-negligible modulation of the X-ray flux but to reach a large RMS am- 
plitude, the blob needs a large azimuthal extent-making it more banana or 
spiral shaped. The reason is that a localized structure does not provide much 
shadowing over a disc rotation period and thus offers only a weak RMS modu- 
lation. In contrast to spirals, such blobs probably cannot account for observed 
LFQPOs in microquasars. 
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4 Expected dependence of RMS amplitude on inclination and en- 
ergy 

Because the modulation comes from shadowing, the disc inclination angle 6 
is important in determining the maximum shadowing and thus the maximum 
RMS amplitude that a given choice of parameters can produce. Motivated by a 
comparison with the LFQPOs of GRS1915+105 (see e.g. the review Mclintock 
& Remillard 2004), we again focus on the parameter choices of simulation #8 
and vary the inclination angle to obtain different values of the RMS amplitude. 
Fig. 5 shows the result. As expected, the more edge-on the view, the higher 
the RMS amplitude. Present observational data are insufficient to definitively 
confirm or contradict the predicted behavior. Such a trend could explain the 
extremely weak LFQPO in Cyg X-l, as the RMS amplitude expected from 
the inferred inclination angle is very small 1 . More objects with a wider range 
of inclination angles are needed. Additional techniques of determining the disc 
inclination angle in microquasars besides using the jet propagation direction 
would also be desirable. In the case of neutron stars, data from eclipsing 
binaries and dippers provides good information on the binary and thus disk 
inclinations (e.g Frank et al., 1987). 

It is important to note that the dependence of the rms amplitude on the 
inclination angle shown in Fig. 5 is displayed only for the case in which all 
other parameters are identical. Therefore the trend shown should be seen 
inside the same "family" of objects-those having similar disk properties. It 
was shown in table 1 that some parameters can have a greater influence on 
the rms amplitude than the inclination. These paramaters could give rise to 
low inclination source with a strong spiral having a higher rms amplitude than 
a high inclination source with a very small spiral 

Using the relation between T and h, we can also study how the RMS amplitude 
behaves as function of energy. Several studies show that LFQPO amplitudes 
rise with energy up to 15-20keV, although the behaviour at higher energies 
seems to depend on the source, for example in GRS 1915+105 (Tomsick & 
Kaaret, 2001). We have not included a corona (high energy component) in 
our simple spectral modeling herein, but we can make a qualitative predic- 
tion of the LFQPO amplitude vs. energy trend by studying the RMS ampli- 
tude as function of our predicted disc flux, supplemented by knowledge of the 
disc/corona flux ratio taken from observation. We find that the RMS ampli- 
tude of the modulation gets bigger at higher energy, his is simply because 
a spiral with r c /rj„ ~ 2 is hotter than the surrounding disc and therefore 
its blackbody contribution to the spectrum peaks at a higher temperature. 



Cyg X-l seems to show a weaker and broader structure than the "usual" LFQPO 
that could be a weak LFQPO in aggrement with this simple model. 
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When a corona is present, this trend would apply at low energies, but we 
expect a critical energy above which the corona strongly dominates and the 
trend reverses. 

We emphasize that in order to predict the exact position of the modification 
in the trend at high energies, we need to take into account not only the disk 
but also the corona self-consistently. This is beyond the scope of the present 
paper which is intended to show the basics of how flux can be modulated by 
non-axisymmetric structures. 



5 Summary and Discussion 

We have presented a simple analytical framework for modeling and interpret- 
ing how non-axisymmetric disk features such as spirals or blobs can create 
modulation of the observed flux from an accretion disk. We connect the disk 
geometry to the flux via the hydrostatic equilibrium assumption. Our focus 
on spirals is motivated by their production in disc simulations of the AEI and 
the efficacy with which such spirals seem to account for the observed LFQPO 
properties of microquasars, see e.g. the review McClintock & Remillard, 2004. 
Our formalism does not depend on the origin of the non-axisymmetric struc- 
tures. Rather, it is used to quantify the influence of a range of generic spirals 
and hotspot geometries. 

Using our formalism (with parameters guided by a fit to the LFQPO proper- 
ties of GRS1915+105) we expect that the RMS amplitude increases with disc 
inclination because the modulation from a spiral or blob results from shad- 
owing of the inner disc. Because the spiral is the hottest part of the disc for 
typical r c ~ 2r in , we also find that the RMS amplitude should increase with 
energy, except above energies where the emission becomes strongly corona 
dominated. More quantitative modeling and more observations are needed. 

If the time scale of a LFQPO in a black holes system is indeed determined 
by an orbit or pattern speed of a non-axisymmetric feature, then at a given 
number of gravitational radii, the LFQPO frequency should scale inversely 
with the black hole mass. If accretion disks around AGN and microquasars 
incur similar instabilities, then a 1 — 10Hz LFQPO in microquasars would 
correspond to a < 10 _6 Hz LFQPO in AGN. Testing this prediction is presently 
difficult because it pushes the present limits of continuous observation times 
for individual AGN. 

Finally note that although we have focused on black hole systems because of 
the LFQPO observations in microquasars, the analytic formalism for modeling 
the influence of spirals and hot spots herein does not depend on the compact- 
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ness of the central object. It can also be applied to discs around neutron stars, 
white dwarfs, or young stellar objects. 

More work is needed to incorporate our basic framework into detailed models 
of individual LFQPO sourceds. 
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Table 1 



For each simulation, all parameters, and the resulting RMS amplitude are given. The 
inclination angle for all cases is 70°, which is the value for GRS1915+105 inferred 
from the orientation of the disk (Mirabel & Rodriguez 1994), assuming the jet is 
perpendicular to the jet. The angular extend of the spiral/blob is given in fraction 
of 2vr. 
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Fig. 1. (a) Schematic top view of of spiral spines (line tracing the peak in height of 
the spiral) for two different choices of a. (b) Schematic slices of the disc height profile 
for an example spiral as a function of radius for h /r = 0.01. Profiles at different 
azimuths are all projected into the same r — z plane. The peak in the left-most 
curve occurs at r c . The width of the curves is determined by the parameter S. The 
downward slope of the line connecting the peaks is determined by (3 (for a fixed a) 
and the maxium height of the peaks is determined by 7 (see text), (c) 3-D close 
up of the inner part of the disc with the spiral wave, the viewing angle is 70° with 
h /r = 0.01, r c = 3, a = 0.07, 7 = 0.3, (3 = 2,5 = 0.2. 
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Fig. 2. Evolution of the RMS as a function of the inclination angle 9 for the 
simulation case #8 of Table 1. 
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